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^\ ' Abstract. Effects of rotational mode coupling on photometric parameters of stellar oscillations are studied. At 

moderate rotation rates, a strong coupling between modes of spherical harmonic degree, £, differing by 2 and 
of the same azimuthal order, m, takes place if the frequencies are close. This is a common situation amongst 
main sequence pulsators. Numerical results for a sequence of f3 Cephei star models are reported for the two- and 
three-mode couplings. 

One consequence of mode coupling is that modes of higher degree should be considered in photometric mode 
identification. Modes with nominal degree £ > 2 acquire substantial £ < 2 components and therefore are more 
likely to reach detectable amplitudes. Coupled mode positions in the amplitude ratio - phase difference diagrams, 
Oh' based on multicolour photometry, become both aspect- and m-dependent. Examples of the mode path in the 

diagram with varying aspect are given. The diagrams remain a useful tool for mode identification in rotating 
s_i ' stars but the tool must be used with care. 
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^ i 1. Introduction Balona & Evers (1999), Cugier & Daszynska (2001), 



Balona et al. (2001) and also by Townsend (2002), who 
applied them to SPB stars. 



Mode identification, that is, determination of the radial 
order and spherical harmonic, is an essential step in aster- 
oseismology. The task is not easy in the case of the oscil- U P to now > the amplitudes and phases, which we - 
lation frequency spectra in (3 Cep and 6 Set stars, which following Cugier et al. (1994) - call photometric nonadi- 
are most often lacking equidistant patterns. The photo- abatic observables, were calculated in the framework of 
metric diagnostic diagrams, i.e., the amplitude ratio vs. linear nonadiabatic theory, ignoring effects of rotation, 
phase difference dependencies in different passbands, are However, amongst p Cep and 8 Set stars slow rotators, 
the most popular tools for mode identification in pulsating for which such an approximation is adequate, are more an 
stars. Following pioneering works (Balona & Stobie 1979, exception than a rule. Here we examine effects of moderate 
Stamford & Watson 1981) these tools have been applied rotation on theoretical diagnostic diagrams. By moderate 
mainly to P Cep and S Set variables. we mean so slow that perturbational treatment of rota- 
Theoretical diagnostic diagrams are based on linear tion is a dequate. Specifically, we rely on the third order 
nonadiabatic calculations of stellar oscillations and on formalism of Soufi et al. (1998). 

models of static plane-parallel atmospheres. In the early The most important effect of rotation in the context 

works an arbitrary parametrization has been used instead of diagnostic diagrams is coupling between close frequency 

of linear nonadiabatic calculations. This approach has modes of spherical harmonic degree , £, differing by 2. The 

been followed even in some recent studies (e.g. Garrido effect was discussed in some detail by Soufi et al. (1998). 

2000). The nonadiabatic calculations were first included The essential formulae are recalled in the next section of 

explicitly by Cugier et al. (1994) and subsequently by this paper. 

Send offprint requests to: J. Daszyhska-Daszkiewicz, e-mail: Numerical results presented later concern one selected 

daszynsk@camk.edu.pl sequence of P Cep models. On a qualitative level the re- 
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suits are applicable to all stars of this type. Our choice of 
(3 Cep stars is motivated not only by the abundance of the 
observational data but also by the fact that we have cred- 
ible results from linear nonadiabatic calculations. This is 
not true for 5 Set stars where there are serious uncertain- 
ties related to the treatment of convection. 

Properties of unstable modes in the selected sequence 
of models are reviewed in Section 3. Also in this section 
we discuss the occurence of near resonances between two 
and three modes as well as certain consequences of mode 
coupling. 

In Section 4 we discuss the visibility in various pass- 
bands of modes described by a single spherical harmonic 
over a wide range of £. Examples of diagnostic diagrams 
for coupled modes are given in Section 5. 

2. Rotational mode coupling 

Chandrasckhar & Lebovitz (1962) were the first to no- 
tice that rotation, even at its slow rate, significantly in- 
fluences mode properties if there is a small frequency dis- 
tance between modes with the same azimuthal order, to, 
and spherical harmonic degree, £, differing by 2. In this 
case, each of the coupled modes must be represented by 
a certain superposition of spherical harmonics of the two 
modes involved. The effect was invoked to explain the non- 
radial character of pulsation in (3 Cep stars. The effect of 
mode coupling has been subsequently discussed in a num- 
ber of papers (e.g. Dziembowski & Goode 1992, Soufi et 
al. 1998). 

In the present paper we assume uniform rotation. 
Although our calculations were carried out to cubic or- 
der in the ratio of rotation to pulsation frequency, Cl/u>, 
for the purposes of this discussion and for the sake of sim- 
plicity, we will write formulae which are accurate only to 
{Q/to) 2 . This is enough to capture the effects we want to 
review here. 

In general, there may be more than two coupled modes. 
Then the displacement eigenfunction of an individual 
mode must be given as a sum 



(i) 



where £o/c stands for the eigenfunctions in the spherical 
model. We choose the normalization in such a way that 
the radial component of the vector £ fc is given by 

Zo, r = RYr^K (2) 

The orthogonality condition is expressed in the form 

d 3 xp£ j ■ £ofc = Sjkh, (3) 



where Ik is called mode inertia. The eigenfuncions and 
eigenvalues for coupled modes are calculated as a solution 
of the matrix equation, 



where E is a unit matrix and a = (ai,a 2 ,...) (see e.g. 
Soufi et al. 1998). The nondiagonal elements, Bjk, are of 
the order of SI 2 and are non-zero for the modes satisfying 
conditions £j = £k ± 2 and mj — nik- With our normal- 
ization these elements are not symmetric. The following 
relation is fulfilled: 



Bkj_ 

Bj k 



1L 
Ik 



(5) 



The diagonal elements are given by 

B kk = (uj ok + mC k nf + 0(fl 2 ), (6) 

where uj^k is the frequency in the spherically symmetric 
model and Ck is the Ledoux constant. The 0(Vl 2 ) term 
arises from quadratic effect of the Coriolis force and from 
the centrifugal distortion. The latter effect dominates for 
nonradial p-modes, as well as for g-modes of low orders. 

Here we limit ourselves to couplings which involve up 
to three modes. The solution in the case of two-mode cou- 
pling is 

2 _ B n + B 22 ± v/(S n - B 22 ) 2 + 4B C 



where Bq 



Bi 2 B 21 , and 
'a-i 



B 



12 



(7) 



(8) 



,G12 / ± B\\ — W ± 

The coupling strength, which is measured by Be, depends 
on mode properties. Coupling between acoustic modes is 
stronger than that involving one or more gravity modes. 
This is so because the effect of the centrifugal distortion 
is only important in the acoustic cavity and it increases 
with the mode frequency. 

For an exact resonance, that is when Bn = B 22 , we 

get 



(9) 



(B - w 2 E)a = 0, 



(4) 



Oi 

, a 2/ ± 

and, regardless of the coupling strength, we have strong 
mode mixing. Then, the relative amplitudes are deter- 
mined by mode inertiae. Therefore, for instance, an I = 2 
mode trapped in the interior thus having a large inertia 
may manifest itself as a radial mode. A close doublet of 
"radial" modes may be detected in this case. 

Now we consider three-mode interaction, adopting 
l 2 = l\ + 2 and £ 3 = £ 2 + 2. An interesting situation 
arises when the closest resonance occurs between modes 
with £ differing by 4. Since we have B\ 3 = B31 = 0, the 
interaction takes place through the l 2 mode. Let us con- 
sider a simple limiting case when B\\ = B 33 . Then, one 
solution is 

lo = B n = B 33 , a 2 = 0, — = -— — . 

The remaining two solutions are the same as in the case of 
two mode coupling, except that Be = B i2 B 2i + B 23 B 32 . 
The highest degree mode, say £ = 4 or 5, which is normally 
undetectable, acquires a large £ = or, respectively, I = 1 
component. This means that some higher degree modes 
may become detectable by means of ground-based pho- 
tometry. 



(10) 
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3. Close frequencies in a sequence of (3 Cephei 
star models 

As an illustration we consider the evolutionary sequence 
of 12 Mq star models in the (3 Cep instability strip. We 
adopted the standard chemical composition (X = 0.7, Z = 

0. 02,. Models were calculated not allowing for convective 
overshooting from the core. In these models we took into 
account the average effect of centrifugal force, assumed 
uniform rotation with an equatorial velocity of 100 km/s 
on the ZAMS, and imposed global angular momentum 
conservation during evolution. The adopted value of rota- 
tional velocity may be regarded typical for (3 Cep stars. 
The coupling strength is proportional to the square of the 
velocity, but the result depends also on the frequency dif- 
ference. Thus, we have no simple scaling of the effect of 
rotational mode coupling. 

Fig. 1 shows the frequency evolution in unstable low- 
degree modes which are designated in accordance with 
the avoided-crossing principle. The physical nature of the 
modes, that is, the relative proportion of the contribution 
from acoustic and gravity propagation zones to mode's 
energy, is reflected in the slope of the cr(logT c ff) function. 
Slow and rapid rises correspond to dominant acoustic and 
gravity mode characters, respectively. The latter implies 
larger mode inertia. 

We have to worry about mode coupling when the fre- 
quency distance between the modes becomes compara- 
ble to the rotational frequency. With our choice of ro- 
tation rate we have dimensionlcss rotational frequency 
cr ro t = fl/^/AirG <p> ~ 0.1. Thus, we can see that in 
many instances the coupling may be significant. 

Examples of coupled mode solutions are given in Table 

1, where only zonal modes (m = 0) are included. The 
values of £' were assigned using again the avoided cross- 
ing principle, that is, the ordering of modes in frequencies 
was kept unchanged. The relevant frequency is vnc, which 
corresponds to B^k ■ In most cases the value of £' identifies 
the dominant component in the surface amplitude. The 
exceptions occur for modes with drastically different iner- 
tiae. 

The cases of exceptionally close two-mode resonances 
are 1, 3, 4, 7, 8, 9 and 11. There are differences between 
these cases, which reflect not only the differences in the 
proximity of the resonance but also in the nature of the 
coupled modes. As we have already discussed in section 

2, the coupling involving modes of gravity character is 
weaker than that involving only modes of the acoustic 
character. Weak coupling is reflected in the low values of 
A^c in cases 1 and 11. However, the mixing of the modes 
is usually strong, except for £' — 1 mode in case 11, which 
remains essentially pure £ = 1. This has to do with the 
disparity of inertiae for the coupled modes. In cases 3 and 4 
both modes have p-mode character and the correction due 
to coupling is larger. Still, the frequency shift of <~ 0.015 
d _1 is significant only for data obtained with observations 
lasting longer than 2 months. In general, we find rather 



small rotation induced frequency changes, even in those 
cases when the coupling is large. 

In some of the listed cases, like 2, 10 and 12, the cou- 
pling appears quite weak. Nonetheless, as we will see later 
in section 5, the diagnostic diagrams may be significantly 
perturbed. 

Among the three-mode couplings, case 9A most closely 
corresponds to the situation described in section 2. With 
the accuracy adopted, the modes £ = 1 and 5 are in the 
exact resonance and the £ = 3 mode is not far away. Also in 
cases 3 A and 7A, the three-mode coupling is important. In 
the latter case we see that £' = 4 mode acquires a large £ = 
component. In case 10A the frequency distance between 
1=1 and 5 is small, but because the distance to £ = 3 is 
relatively large, the coupling is essentially inconsequential. 
Also in case 6A, the coupling to the third mode is without 
significant consequences. Here the reason is the gravity 
character of £ — 4 mode. 

Coming to tesseral modes (to ^ 0), we first note that 
the coupling occurs if m = ±1 and, for even £ modes, also 
if m = ±2. In general, coupling conditions between pro- 
grade and retrograde mode pairs are not the same. The 
asymmetry is due to the difference in the Ledoux constant 
for the coupled modes. The only systematic difference rel- 
ative to the to = modes is a somewhat weaker coupling 
due to smaller value of the angular integrals entering ex- 
pressions for Bjk- 

The near degeneracy of rotationally coupled modes oc- 
curs also in wide ranges of 5 Set star models. For high 
order p-modes the effect is systematic, as implied by the 
asymptotic relation cj n j+2 ~ w n -i,e- Since the coupling 
strength increases with frequency, in the high order p- 
mode pulsators the effect begins to be significant at lower 
rotation rates than in (3 Cep stars. Frequency changes 
caused by rotational mode coupling in 5 Set stars and in 
solar type pulsators have been discussed by Goupil et al. 
(2000) and by Dzicmbowski & Goupil (1998), respectively. 

4. Photometric diagnostic diagrams for single 
spherical harmonic modes 

Our approach follows that of Cugier et al. (1994), ex- 
cept for correcting an error in the expression for per- 
turbed gravity, which fortunately has negligible conse- 
quences (Cugier & Daszyhska, 2001). Here we outline the 
main steps. In our construction of theoretical diagnostic 
diagrams we rely on the linear nonadiabatic description of 
oscillations. Normalization of eigenfunctions is fixed in Eq. 
2. With that, the local displacement of the photosphere in 
the stellar reference frame may be written as follows 

Sr(R, 9, cp) = eRRe{Y e m e - iojt }, (11) 

where e is a small quantity, which is determined by data. 
This equation implies that the phase zero corresponds to 
the maximum of displacement. The use of —icot in the ex- 
ponential time dependence, which is now the most com- 
mon choice, instead of +itot used by Cugier et al. (1994), 
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Fig. 1. Unstable low degree modes in a sequence of j3 Cephei star models with mass of 12 M©. The models cover only 
the expansion phase of the main sequence evolution. In the left panel note that the £ = 0,p\ and £ = 2, g\ modes have 
close frequencies in the extended range of logT e g. At lower logT e ff the I = 2,g 2 mode comes to a near resonance with 
the same radial mode. We see also the proximity of the nonradial I = 2 and 4 modes at a w 1.6 — 1.8. In the right 



panel, note a near resonance between the I 
three modes occurs only occasionally. 



1,01 mode and, first, the I = 3,02 and, later, the g$. The proximity of 



results in the opposite signs for the phase differences be- 
tween calculated observables. For a comparison of calcu- 
lated and measured values it is important to check the 
sign convention adopted in the data analysis. With the 
sign convention adopted in the present work, the phase 
difference 4> a — (j>b > means that the maximum of a oc- 
curs after the maximum of b. 

The corresponding changes of the effective tempera- 
ture and gravity during the pulsation cycle are given by 



where 



and 



^ffeff 
9cS 



3uj 2 \ 6r 
4ttG <p> J ~R' 



(12) 



(13) 



respectively, where / is the complex quantity determined 
from linear nonadiabatic calculations, < p > is the mean 
density of the star and G is the gravitational constant. 

Further, we assume static plane-parallel atmosphere, 
and we write the monochromatic flux variation caused by 
the oscillation mode with £ degree in the following form: 



A^ A 

A 



sY e m ^0)b^Re{[Dl e + D 2A + DUe^}, (14) 



1 d\og(T x \b x t \ 
4 J aiogr off 



D 2 ,t = (2 + £)(!-■ 



D 



AttG <p> 



dlog(F x \b}\) 



and 



/' 

Jo 



h°x{p)p,Pt{p)d^, 



(15a) 
(156) 
(15c) 

(16) 



In this formula, the D x e term describes the temperature 
effects, whereas the influence of the gravity changes is 
contained in the D% e term. Both include the perturba- 
tion of the limb-darkening. In our calculation we relied 
on Claret's (2000) analytical fit to Kurucz's (1998) tabu- 
lar data. The ^-dependence of Di and D 3 arises from the 
nonlincarity of the adopted limb-darkening law. The D 2y i 
term stands for the geometrical effects. The effect of ori- 
entation with respect to the observer is described by the 
spherical harmonic Y™(i, 0). 

In Fig. 2 we show the dependence of the disc averag- 
ing factor, bp, on mode degree, £, for different passbands. 
For comparison, the bolomctric values are also shown. 
These are similar to those calculated by Dziembowski 
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Table 1. Selected cases of rotational mode coupling. Co is the dimensionless oscillation frequency in the spherical 
models, v a — is the corresponding cyclic frequency in d _1 ; u^c takes into account effects of rotation except for 

mode coupling; Ave is a frequency correction due to coupling, whereas Avt is a total frequency change relative to 
the spherical model; au are amplitudes of the involved eigenmode components. 
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coupling 


ffo 
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modes 
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CS2 


03 


1 


4.374 


l\ = 2, g 2 


1.611 


3.836 


3.849 


0.002 


0.014 


£' 


= 2 


0.88 


-0.47 








(■2 = 4, g 4 


1.610 


3.834 


3.846 


-0.002 


0.011 


£' 


= 4 


0.62 


0.78 




2 


4.406 


£i = 0, pi 


1.875 


6.438 


6.457 


0.002 


0.022 


£' 


= 


1.00 


-0.07 








£2 = 2, 9l 


1.801 


6.182 


6.215 


-0.002 


0.031 


£' 


= 2 


0.14 


0.99 




3 


4.388 


£1 = 0, pi 


1.913 


5.265 


5.282 


0.013 


0.030 


£' 


= 


0.91 


-0.42 








£2 = 2, 5 i 


1.894 


5.214 


5.248 


-0.013 


0.021 


£' 


= 2 


0.51 


0.86 




3A 


4.388 


£1 = 0, pi 


1.913 


5.265 


5.282 


0.011 


0.028 


£' 


= 


0.92 


0.36 


-0.15 






£2 = 2, gi 


1.894 


5.214 


5.248 


-0.014 


0.020 


£' 


= 2 


0.50 


-0.86 


0.05 






£3 = 4, 52 


1.921 


5.288 


5.303 


0.003 


0.019 


£' 


= 4 


0.59 


0.51 


0.63 


4 


4.374 


£1 = 0, pi 


1.932 


4.600 


4.615 


-0.015 


0.000 


£' 


= 


0.85 


0.53 








£2 = 2, gi 


1.934 


4.605 


4.636 


0.015 


0.046 


£' 


= 2 


0.58 


-0.82 




5 


4.392 


£1 = 0, P 2 


2.462 


7.096 


7.110 


0.008 


0.022 


£' 


= 


0.99 


-0.16 








£2 = 2, po 


2.385 


6.877 


6.913 


-0.008 


0.028 


£' 


= 2 


0.24 


0.97 




6 


4.376 


£1 = 0, p 2 


2.496 


6.064 


6.076 


0.015 


0.027 


£' 


= 


0.95 


-0.30 








£2 = 2, po 


2.445 


5.939 


5.976 


-0.015 


0.022 


£' 


= 2 


0.39 


0.92 




6A 


4.376 


£1 = 0, p 2 


2.496 


6.064 


6.076 


0.015 


0.027 


£' 


= 


0.95 


0.30 


0.00 






£2 = 2, po 


2.445 


5.939 


5.976 


-0.015 


0.022 


£' 


= 2 


0.39 


-0.92 


0.00 






£3 = 4, gi 


2.491 


6.050 


6.058 


0.000 


0.009 


£' 


= 4 


0.17 


-0.06 


-0.98 


7 


4.354 


£1 = 0, p 2 


2.508 


4.875 


4.886 


-0.019 


-0.009 


£' 


= 


0.91 


0.41 








£2 = 2, po 


2.521 


4.899 


4.925 


0.019 


0.046 


£' 


= 2 


0.59 


-0.81 




7A 


4.354 


£1 = 0, p 2 


2.508 


4.875 


4.886 


-0.020 


-0.009 


£' 


= 


0.91 


-0.42 


0.00 






£2 = 2, po 


2.521 


4.899 


4.925 


0.020 


0.046 


£' 


= 2 


0.59 


0.81 


0.01 






£3 = 4, 52 


2.518 


4.895 


4.902 


-0.000 


0.007 


£' 


= 4 


0.71 


0.27 


-0.64 


8 


4.370 


£1 = 1, 5i 


1.981 


4.532 


4.566 


-0.023 


0.011 


£' 


= 1 


0.73 


0.68 








£2 = 3, 52 


1.983 


4.537 


4.571 


0.023 


0.057 


£' 


= 3 


0.67 


-0.75 




8A 


4.370 


4 = 1, gi 


1.981 


4.532 


4.566 


-0.018 


0.016 


£' 


= 1 


0.79 


-0.60 


-0.12 






£2 = 3, 52 


1.983 


4.537 


4.571 


0.027 


0.061 


£' 


= 3 


0.61 


0.79 


0.10 






£3 = 5, 54 


1.940 


4.439 


4.456 


-0.009 


0.009 


£' 


= 5 


0.09 


-0.42 


0.90 


9 


4.360 


*i = 1, 9i 


2.006 


4.144 


4.180 


0.022 


0.058 


£' 


= 1 


0.75 


-0.66 








£2 = 3, 52 


2.008 


4.148 


4.178 


-0.022 


0.007 


£' 


= 3 


0.72 


0.70 




9A 


4.360 


ft = 1, 9i 


2.006 


4.144 


4.180 


0.000 


0.036 


£' 


= 1 


0.87 


0.00 


-0.49 






£2 = 3, 52 


2.008 


4.148 


4.178 


-0.044 


-0.014 


£' 


= 3 


0.38 


-0.71 


0.59 






£3 = 5, 54 


2.006 


4.143 


4.180 


0.044 


0.081 


£' 


= 5 


0.39 


0.69 


0.61 


10 


4.406 


ft = 1, pi 


2.021 


6.938 


6.988 


0.004 


0.055 


£' 


= 1 


1.00 


-0.08 








£2 = 3, gi 


1.927 


6.616 


6.663 


-0.004 


0.043 


£' 


= 3 


0.15 


0.99 




10A 


4.406 


£1 = 1, Pi 


2.021 


6.938 


6.988 


0.002 


0.053 


£' 


= 1 


0.97 


0.05 


-0.25 






£2 = 3, 51 


1.927 


6.616 


6.663 


-0.013 


0.034 


£' 


= 3 


0.14 


-0.96 


0.23 






£3 = 5, 51 


2.022 


6.940 


7.007 


0.011 


0.077 


£' 


= 5 


0.21 


0.12 


0.97 


11 


4.372 


£1 = 1, Pi 


2.230 


5.205 


5.224 


0.000 


0.019 


£' 


= 1 


1.00 


0.05 








£2 = 3, gi 


2.231 


5.208 


5.218 


-0.000 


0.009 


£' 


= 3 


0.66 


-0.75 




12 


4.380 


£l = 1, P2 


2.690 


6.795 


6.858 


0.009 


0.072 


£' 


= 1 


0.99 


-0.13 








£2 = 3, po 


2.584 


6.528 


6.572 


-0.009 


0.035 


£' 


= 3 


0.21 


0.98 




12A 


4.380 


£l = 1, P2 


2.690 


6.795 


6.858 


0.009 


0.072 


£' 


= 1 


0.99 


0.13 


0.00 






£2 = 3, po 


2.584 


6.528 


6.572 


-0.009 
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(1977) with the Eddington's limb darkening law. In con- 
trast, large differences are seen for individual passbands. 
Stronger averaging of odd-£ modes gives rise to a pro- 
nounced oscillatory dependence on £, which is superim- 
posed on the decreasing trend. The amplitude is larger 
than for the bolometric flux due to the natter limb- 



darkening law. Note that in the case of uniform brightness 
disc, the odd-£ modes beginning with I — 3 are completely 
averaged out. 

The decreasing trend does not properly reflect the I de- 
pendence of the expected amplitudes. This would be true 
only if the temperature term dominates. For modes of low 
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Fig. 2. The absolute values of the disc averaging fac- 
tor plotted for bolometric and for twelve passbands as 
a function of £. The passbands include u,v,b,y and 
U, B, V, R, I, J, H, K filters. 

radial orders the geometrical term D 2 becomes dominant 
around I = 6. Then owing to the (£ + 2)(£ — 1) factor the 
decline with £ is slower than implied by Fig. 2. 

For calculating diagnostic diagrams it is convenient to 
use the amplitude of the monochromatic flux variation in 
the following complex form 



A x (i) = e y/"(i,0)#(£>k + D 2 j + D* t ), 



(17) 



which is equivalent to Eq.14. From this expression we can 
directly obtain the amplitude ratio and phase difference 
in any selected pair of passbands. Our choice is u and y 
Stromgren filters. For modes described by a single Y™, 
these photometric observables are independent of the as- 
pect angle and azimuthal order, m. 

The diagnostic diagram for unstable modes in our 
model sequence is presented in Fig. 3. Modes up to £ = 6 
are included. We can see that, in contrast to the radial 
modes, the nonradial modes are crowded in small areas. 
Two aspects of this plot call for an explanation. Firstly, 
why it is that the amplitude ratios and phase differences 
of radial modes are particularly dependent on frequency 
and model temperature? Secondly, what is the origin of 
the grouping of nonradial modes? 

The main reason for this distinctive behavior of ra- 
dial modes is a comparable but opposite contribution of 
the temperature and the geometrical terms to luminosity 
variation. The sum of the real part of D\ and D 2 varies 
significantly with temperature and it is strongly frequency 
dependent. For £ = 1, we have D 2 = and the only reason 
for the spread of the mode positions is a weak dependence 
of b\ on the passbands. The role of the geometrical term 
for £ = 2 modes is comparable to that of D*. However, 
the two terms add, which results in much smaller spread 
of their positions. The near overlap of I = 1 and 3 modes 
is coincidencial. At still higher £, when geometrical effect 
dominates, the even and odd degree modes are gathered 
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Fig. 3. Theoretical diagnostic diagram with the 
Stromgren passbands u and y for unstable low degree 
modes in the (3 Cep star models with M — 12M Q . In the 
sequences for the I = modes the most negative values of 
the phase difference correspond to the coolest models. 



in two small separate domains, determined solely by the 
b\ factor (see Fig. 2). The phase difference is always close 
to zero owing to the small role of the only complex term, 
D\. 



5. Photometric diagrams for coupled modes 

The complex amplitude of the monochromatic flux varia- 
tion for a coupled mode may be expressed in the form 



(18) 



where are solutions of Eq. 4 and A\(i) is given in Eq.17. 
The quantities at describe contributions of the £fc-modes 
to the coupled mode £'. Selected values of are given 
in Table 1. The quantity / appearing in Eq. 17 is only 
weakly dependent on £. Its relatively strong dependence 
on frequency is irrelevant in view of very small frequency 
spread of the mode involved. In contrast to the ampli- 
tude ratio and phase difference for the pure modes, these 
quantities now depend on the inclination angle, i, and the 
azimuthal order, m. We will show here few representative 
examples of the photometric diagnostic diagrams in the 
case of two and three mode coupling and we will discuss 
the consequence for the mode identification. 

Naively, one could expect that the coupled mode ap- 
pears in the diagnostic diagrams between positions cor- 
responding to its components associated with individual 
spherical harmonics. We will see that this is not true. 

As our first illustrative example, we chose the case 
when the £ = and 2 modes are in close resonance (case 4 
in Table 1). In Fig. 4 we show how positions of the coupled 
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Fig. 4. The diagnostic diagrams with Stromgrcn pho- 
tometry passbands u and y, showing positions for cou- 
pled I = 0,pi and £ = 2,gi modes (case 4 in Table 1) in 
the model at logT cff = 4.3741 (see Fig. 1). This model 
corresponds to the smallest frequency distance between 
these two modes. The upper panel refers to the solution 
dominated by the I — component, while the lower one 
refers to the solution dominated by £ = 2 component. 
Arrows correspond to observations from the polar direc- 
tion. Spacing between consecutive dots is 0.02 in cosi. 
Note that the density of dots reflects probability of ob- 
serving the modes in various parts of that diagram. Circles 
indicate positions of pure I = 0, 1, 2 modes. 



modes move in the diagnostic diagram with varying as- 
pects. A comparison of plots in the upper and lower pan- 
els reveals that the movement is strongly mode-dependent. 
Let us focus first on the case of mode dominated by the 
radial component, shown in the upper panel. Due to close 
resonance, contamination with the £ = 2 component is 
significant. The amplitudes of the £ = and 2 compo- 
nents are 0.85 and 0.53, respectively. An observer from 
the polar direction will identify this mode as I — 1. This 
is so, because the main role of the £ = 2 component is a 
cancellation of the geometrical term in the expression for 
the flux variation (see. Eq. 14). When the observer moves 
away toward the equatorial direction, he will see the mode 
as a pure £ = 0, at the angle corresponding to the node of 
F 2 °(i, 0), and subsequently he will get mode characteristics 
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Fig. 5. Complex amplitudes in Stromgren passbands u, y 
for modes used in Fig. 4. The three pairs of vectors shown 
in the upper panel correspond to cosi = 1.0, 0.32 and 0.0. 
The four pairs of vectors in the lower panel correspond to 
cos i = 0.0, 0.66, 0.72 and 1.0. The arrows mark amplitudes 
seen from the polar direction. 




Fig. 6. The same as in Fig. 4 but the coupling between 
modes £ = and 2 is much weaker (case 2 in Table 1). 



which do not correspond to any of single-^ modes. Note 
that such an appearance has a significant probability, as 
measured by the density of dots. 

Fig. 5 helps us to understand the cause of the described 
behaviour. Here we show complex, arbitrarily normalized, 
amplitudes in the two selected passbands. To make the 
amplitudes realistic one should divide them by a factor 
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between 10 to 100. We see an increasing phase difference 
between two passbands when we move from polar to equa- 
torial direction. This is predominantly due to a decrease 
of the real part of the complex amplitude in y passband. 
Why does the argument for the y passband change much 
more rapidly with aspect than that for u passband? This 
is due to a smaller absolute value of D\ , hence the greater 
role of D 2 - Only in this colour does Re(.4) pass through 
zero. This happens at cos-i = 0.32. Note, however, that 
when the phase difference is large, the pulsation ampli- 
tude is rather small. So that the situation may not be 
easily observable. 

The very different pattern of the £' — 2 mode be- 
haviour, shown in lower panel of Fig. 4, is due to the fact 
that now the real part of the amplitude in both colours 
changes sign. This is caused by the greater contribution 
from the £ = 2 component. Here the equatorial position of 
the mode is not so distant from the polar one. Surprisingly, 
the mode is often seen as either £ = 1 or £ > 3. 

Fig. 6 illustrates case 2, when as we have noted the 
coupling is rather weak. Whereas for the I' = root the 
effect of the coupling is marginal, for I' = 2 it is quite 
significant. In spite of the fact that the mode is dominated 
by £ — 2 component, it has a fairly high chance to be 
identified as £ = 1 or £ > 3. 

The cases of I = 1 and 3 coupling never lead to a large 
departure of the I' = 1 modes from the £ = 1 position. 
Since mode £' = 3 acquires a substantial £ — 1 compo- 
nent, it is mostly seen through this component and hence 
would be identified as I = 1. It is interesting that many 
different modes in rotating stars may be identified as £ = 1 
by means of the diagnostic diagram method, if effects of 
rotation are ignored. 

Fig. 7 illustrates an interesting case (case 3A in Table 
1) of the three- mode coupling. The I 1 = mode is most 
often seen close to the £ = position, but there is also a 
fair probability to see this mode close the £ = 1 position. 
The path for the £' = 2 mode is very similar to that shown 
in the lower panel of Fig. 5. The £' = 4 mode is most 
often seen between £ = and £ = 1 positions. However, 
seen from the equatorial direction it is quite far from any 
position of a single £ mode. 

For other modes of the same multiplets the coupling 
effects are less spectacular. In all the cases from m = —2 
to m = +2, the £' = 2 modes hardly move from the £ = 2 
position. As expected, the effect for the £' = 4 mode is 
more significant. The mode stays predominantly between 
the £ = 1 and I — 4 positions. 

Fig. 8 shows how mode-coupling of the multiplets 
would be identified if rotation is neglected. In this fig- 
ure we plotted the observed amplitudes assuming that the 
intrinsic mode amplitudes are the same. Note that the in- 
dicated amplitudes should be divided by a factor between 
10 and 100 to correspond to reality Identifications were 
based on the diagram in Fig. 3. 
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Fig. 7. An example of the diagnostic diagrams for the 
three mode coupling (case 3A). 



6. Conclusions 

We have seen that even at moderate rotation, mode cou- 
pling leads to complications in the diagnostic diagrams 
used for photometric mode identification. Modes of de- 
grees I > 2, which are often ignored in identification of 
peaks in oscillation spectra, may acquire substantial low I 
components, and are more easily detected. Unlike modes 
described by a single spherical harmonic, the positions of 
coupled modes in the diagnostic diagrams depend on the 
aspect and on the azimuthal order. The positions may 
be quite confusing, for instance, a mode composed of the 
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Fig. 8. Multiplcts of coupled £' = 0,2,4 modes for case 
3A. Amplitudes in the u band were calculated assuming 
common normalization, e = 1 in Eq.17. Numbers above 
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effects of rotation. Lack of numbers means that mode can 
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£ = and £ = 2 components may appear for a range of 
aspect angles at the £ = 1 position. 

All that is not good news for the photometric mode 
identification procedure. The problems are not confined 
to a few cases, but occur at the typical rotation rates 
encountered in (3 Cep and 5 Set stars. Close frequencies 
of rotationally coupled modes occur over wide ranges of 
the instability strips. The implication is that we must be 
careful in using the diagnostic diagrams for inferring the 
£ values. The diagrams remain useful. After all, they do 
provide observational constrains on stars and their oscil- 
lations. However, mode identification may be done only 
simultaneously with determination of stellar parameters 
and inclination of rotational axis. 
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